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Abstract 

The  engineering  of  advanced  optoelectronic  integrated  circuits  implies  the  stringent  control  of 
thickness  and  composition.  These  demands  led  to  the  development  of  surface-sensitive  real-time 
optical  sensors  that  are  able  to  move  the  control  point  close  to  the  point  where  the  growth  occurs, 
which  in  a  chemical  beam  epitaxy  process  is  the  surface  reaction  layer  (SRL),  built  up  of 
physisorbed  and  chemisorbed  precursor  fragments  between  the  ambient  and  film  interface.  In  this 
contents,  we  explored  the  application  of  p-polarized  reflectance  spectroscopy  (PRS)  for  real-time 
monitoring  and  control  of  pulsed  chemical  beam  epitaxy  (PCBE)  during  low  temperature  growth  of 
epitaxial  Ga1.xInxP  heterostructures  on  Si(001)  substrates.  A  reduced  order  surface  kinetics 
(ROSK)  model  has  been  developed  to  describe  the  decomposition  and  growth  kinetics  of  the 
involved  organometallic  precursors  and  their  incorporation  in  the  film  deposition.  We  demonstrate 
the  linkage  of  the  PRS  response  towards  surface  reaction  chemistry,  composition,  film  growth 
rate,  and  film  properties.  Mathematical  control  algorithms  are  applied  that  link  the  PR  signals  to  the 
growth  process  control  parameters  to  control  composition  and  growth  rate  of  epitaxial  Ga,_xInxP 
heterostructures. 
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1.  Introduction 


Applying  optical  probe  techniques  to  real-time  characterization  of  thin  film  growth  inherits  the 
challenge  of  relating  surface  chemistry  processes,  that  drive  the  growth  process,  to  growth/film 
properties,  such  as  composition,  instantaneous  growth  rate  or  structural  layer  quality.  The  need  is 
especially  acute  for  chemical  deposition  methods,  where  the  surface  plays  a  major  role  in  the 
growth  process  and  the  addition  of  small  amounts  of  new  reactants  (e.g.,  dopants)  can  severely 
modify  growth  chemistries1’2.  The  limited  knowledge  in  these  areas  slowed  the  progress  in 
understanding  and  controlling  thin  film  growth.  To  improve  the  understanding  of  the  driving 
mechanisms  of  growth  processes,  non-intrusive  real-time  techniques  have  been  developed, 
focusing  on  the  monitoring  of  surface  processes  by  reflection  high  energy  electron  diffraction 
(RHEED)3,  reflectance  difference  spectroscopy  (RDS)4,  surface  photo  absorption  (SPA)5-7  and  p- 
polarized  reflectance  spectroscopy  (PRS)8-10. 

Presently,  the  only  two  techniques  that  combine  the  advantage  of  high  surface  sensitivity  with 
bulk  film  properties  characterization  are  (a)  an  integrated  spectral  ellipsometry  (SE)  /  RDS 
spectrometer  developed  by  D.  Aspnes  et  al.11  and  (b)  PRS8-10.  Both  techniques  aim  to  integrate 
the  optical  response  to  surface  processes  with  the  optical  response  to  bulk  properties  to  monitor 
and  control  the  deposition  process  with  sub-monolayer  resolution. 

This  contribution  describes  the  recent  developments  utilizing  p-polarized  reflectance 
spectroscopy  (PRS)  for  closed-loop  deposition  control  during  pulsed  chemical  beam  epitaxy 
(PCBE)  of  III-V  heteroepitaxial  growth.  The  demonstrated  high  sensitivity  of  PRS  towards 
surface  reactions  processes  in  the  context  of  real-time  monitoring  of  PCBE  has  opened  new 
possibilities  for  characterization  and  control  of  thin  film  deposition  processes.  During 
heteroepitaxial  GaP/Ga^n^F  growth  on  Si  under  PCBE  conditions  the  surface  is  periodically 
exposed  to  metalorganic  precursors,  which  causes  a  periodic  in  composition  and  thickness  altered 
surface  reaction  layer  (SRL).  The  control  of  a  growth  process  using  the  optical  signature  from  the 
SRL  that  feeds  the  underlying  growth  requires  detailed  instantaneous  simulation  and  prediction  of 
the  surface  chemistry  and  its  link  to  the  optical  properties  of  the  outer  most  layer  in  a  multilayer 
medium.  A  reduced  order  surface  kinetics  (ROSK)  model  has  been  developed  that  describes  the 
growth  process  with  a  mathematically  reduced  number  of  surface  reactions  equations  using 
heteroepitaxial  GaxIni_xP  growth  as  an  example.  The  dynamic  in  the  molar  concentrations  of 
surface  constituents  evolution  gives  information  on  SRL  thickness,  its  optical  response  in  a  four 
media  layer  approximation,  the  instantaneous  growth  rate,  and  the  composition  of  the  growing 
film.  For  real-time  closed-loop  deposition  control  a  virtual  substrate  approach  was  used,  an 
approach  recently  introduced  by  D.E.  Aspnes  for  product-driven  deposition  control12. 

2.  Experimental  Setup  and  Results 

For  monitoring  both  the  bulk  and  surface  properties  during  heteroepitaxial  GaxIni_xP  growth  on 
Si,  p-polarized  reflectance  spectroscopy  (PRS)  has  been  integrated  in  a  pulsed  chemical  beam 
epitaxy  (PCBE)  system  as  schematically  shown  in  Fig.  1(b).  In  PCBE,  the  surface  of  the 
substrate  is  exposed  to  pulsed  ballistic  beams  of  TBP  [(C4H9)PH2]  and  TEG  [Ga(C2H5)3]  and 
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TMI  [In(CH3)3]  at  typically  350  -  450°C  to  accomplish  nucleation  and  overgrowth  of  the  silicon  by 
an  epitaxial  GaP  /  Gai_xInxP  film.  For  PRS  and  laser  light  scattering  (LLS)  we  employed  p- 
polarized  light  beams  (A,i=632.8  nm  and  7.2=700  nm  )  and  Glan-Thompson  prisms,  as  illustrated 
in  Fig.  1(a).  The  beams  impinge  on  the  substrate  at  two  angles  of  incidence  PR70  ((p=71.5°)  and 
PR75  (tp=75.2°).  Further  details  on  the  experimental  conditions  are  given  in  previous 
publications8-10’13-23. 


Figure  1:  (a)  Schematic  setup  of  growth  monitoring  by  PRS,  LLS,  and  QMS;  and  (b)  Setup  of 
PCBE  system  for  III-V  compound  semiconductor  growth 


Figure  2  shows  the  evolution  of  the  PR  signals  during  growth  of  Ga|_xInxP/GaP  on  Si(001)  at 
420°C,  recorded  for  PR70  and  PR75  at  A.=650  nm  ±5nm  and  at  A.=632.8  nm,  respectively.  The 

growth  process  is  composed  in  4  sections: 

(a)  substrate  and  surface  preconditioning; 

(b)  deposition  of  a  GaP  buffer  layer  lattice-matched  to  the  substrate; 

(c)  growth  of  a  graded  Ga^I^P  layer,  which  is  shown  in  more  detail  in  Fig. 3;  and 

(d)  growth  of  a  GaP  cap  layer. 

During  the  pre-conditioning  period,  the  PR  signals  change  according  to  the  temperature 
dependency  of  the  substrate.  The  signals  are  used  to  verify  independent  temperature  measurements 
and  to  cahbrate  the  actual  surface  temperature.  After  initiating  growth,  minima  and  maxima  are 
observed  in  the  time  evolution  of  the  PR  signals  due  to  the  interference  phenomena  as  the  film 
grows.  Please  note  that  both  signals  are  phase  shifted,  which  is  due  to  the  fact  that  one  angle  of 
incidence  (PR75)  is  above  -  and  the  other  (PR70)  below  -  the  pseudo-Brewster  angle  of  the 
growing  film  material.  This  angle  is  defined  for  an  infinitive  thick  layer  only  and  is  in  the  range 
between  71  and  74  deg  for  Gaj.xInxP  monitored  at  632.8  nm.  Superimposed  on  the  interference 
oscillations  of  the  reflected  intensity  is  a  fine  structure  that  is  strongly  correlated  to  the  time 
sequence  of  the  supply  of  precursors  employed  during  the  steady-state  growth  conditions. 
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Figure  2:  Growth  monitored  by  PRS  during  heteroepitaxial  Gaj.xInxP/GaP  on  Si(001). 


Figure  3: 

Growth  of  a  graded  Ga^_ 
xInxP  layer  (see  graded  Ga  |  _ 
xInxP  region  marked  in  Fig. 
2): 

(a)  Estimated  composition, 
determined  via  ex-situ 
x-ray  diffraction  analysis. 

(b)  and  (c):  evolution  of  the 
PR  signals  .  The  insets  show 
the  fine  structure  response  at 
two  different  positions  with 
different  TMFTEG  flow  ratios 
and  different  PR  responses  to 
it.  The  ROSK  model 
introduced  below  will 
demonstrate  how  the  optical 
PR  response  is  linked  to 
composition  and  growth  rate 
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3.  Reduced  order  surface  kinetics  (ROSK)  model  for  Gai-xInxP  deposition 

The  reduced  order  surface  kinetics  (ROSK)  model  for  the  binary  compound  semiconductor  GaP 
growth  from  trie  thy  lgalli  um  and  tertiary-butylphosphine  TBP  [(C4H9)PH2]  has  been  discussed 
previously23.  At  this,  we  extend  this  model  by  adding  the  decomposition  mechanism  for 
trimethylindium  (TMI)  for  the  description  of  heteroepitaxial  GaxIri|_xP  growth.  The  TMI 
defragmentation  and  sufficient  retention  of  fragments  on  the  surface  occurs  within  a  limited  process 
window  in  the  temperature  range  700K  <  T  <  850K  for  Si(100)24.  The  kinetics  of  TMI  pyrolysis 
for  the  growth  of  GaxInj_xP  utilizing  trimethylindium,  triethylgallium  and  tertiary-butylphosphine 

as  source  vapors  has  been  discussed  in  detail  elsewhere20’25.  Its  progression  can  be  summarized  in 
three  consecutive  steps: 


In(CH3)3  In(CH3)2*  +  CH3‘, 

(1) 

In(CH3)2  ->  HnCH3  +  CH3‘, 

(2) 

IInCH3  lln*  +  CH3*, 

(3) 

where  the  vertical  dashes  and  superscript  dots  denote  lone  electron  pairs  and  single  valence 
electrons,  respectively. 

For  GaP  growth  on  Si(100),  we  have  shown  that  the  decomposition  of  TBP  is  fast  and 
elimination  of  ethyl  radicals  from  the  TEG  fragments  represents  the  rate  limiting  step26.  Depending 
on  the  delay  between  the  TEG  and  TBP  source  vapor  pulses,  carry-over  of  TEG  fragments  from 
one  precursor  pulse  cycle  to  the  next  may  occur,  which  establishes  in  steady-state  a  surface 
reaction  layer  (SRL),  consisting  of  a  mixture  of  reactants  and  products  of  the  chemical  reactions 
that  drive  the  epitaxial  growth  process.  In  a  realistic  model,  the  SRL  represents  a  multicomponent 
mixed  phase  with  a  variety  of  radical  reactions  that  have  to  be  added  to  the  above  reactions  (1)  -  (3) 
and  to  the  reaction  products  from  the  TEG  and  TBP  decomposition26.  The  thickness  and 
composition  of  the  SRL  depends  on  the  relative  heights  and  widths  of  the  employed  TMI,  TEG 
and  TBP  source  vapor  pulses  and  their  repetition  rate.  We  note  that  some  of  the  intermediate 
fragments  of  the  source  vapor  molecules  in  the  SRL  that  feed  III-V  CBE  may  carry  permanent 
dipole  moments,  which  are  likely  to  contribute  to  the  stabilization  of  the  SRL.  In  view  of 
intermolecular  interactions,  deviations  of  the  SRL  from  ideal  behavior  can  be  expected.  However, 
the  objective  here  is  to  relate  the  measured  PR  signals  to  the  dynamic  of  the  decomposition 
processes  of  the  SRL  constituents,  and  to  relate  the  kinetics  of  growth  and  real-time  modeling  for 
closed  loop  process  control. 

The  reduced  order  kinetic  model  for  the  compound  semiconductor  Ga|_xInxP,  summarizes  all 
chemical  reactions  in  one  dominant  bimolecular  reaction  for  the  TBP  pyrolysis  (first  precursor, 
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PCI),  two  dominant  bimolecular  reactions  for  the  TEG  decomposition  (second  precursor,  PC2), 
and  two  dominant  reactions  for  the  TMI  decomposition  process.  All  precursors  are  supplied 
sequentially  separated  by  pauses  as  shown  schematically  in  Fig.  4 


0  0.8  t]  1 1+0.3  3.0  3.8  t2  t2+0.3  6.0 

time  (sec) 


Figure  4:  Schematic  representation  of  a  precursor  cycle  sequence  used  for  the  growth  of  the 
ternary  compound  semiconductor  Gai_xInxP  grown  via  the  organo-metallic  precursors 
TBP,  TEG  and  TMI.  The  start  positions  for  TEG  and  TMI,  indicated  by  tj  and  t2  were 
1.5  sec  and  4.5  sec  for  the  experiment  shown  in  Fig.  2. 


With  the  above  mentioned  assumptions,  the  differential  rate  equations  for  the  molar 
concentrations  nj  of  SRL  constituents  in  the  Gaj.xInxP  system  can  be  written  as 


—*4(0=  nXBP  -  nj(t)  -  jf4  n3(t)  nj(t)  -  jf7  n6(t)  nj(t), 
d  t 


—  n2(t)  -  nTEG  -  £2  n2(t)’ 
d  t 


— n3(t)=  i-2  n2(t)  -  £3  n3 (t)  -  §f4  n3(t)  nj(t), 
d  t 


—  n5(t)  -  nTMI  -  n5(t)  , 

d  t 


—  n6(t)=  sf5n5(t)  -  4n6(t)  -  8^7  n6(t)  nj(t), 
d  t 

with  the  two  incorporation  reactions 


—  n4(t)=  #4  n3(t)  nt(t). 
d  t 


(4) 

(5) 

(6) 

(7) 

(8) 

(9) 


6 


y-n7(t)  =  in  n6(t)  n^t). 
d  t 


(10) 


for  GaP  and  InP,  respectively. 


Equation  (4)  through  (6)  describes  the  reduced  order  TBP  and  TEG  pyrolysis  as  described 
earlier23,  n-pgp  and  iijjtq  denotes  periodic  supply  functions  expressed  in  terms  of  the  molar 
concentration  of  TBP  and  TEG  reaching  the  surface.  To  this,  we  add  Equation  (7)  and  (8)  that 
describe  the  parameterized  reduced  order  TMI  pyrolysis.  This  pyrolysis  is  assumed  to  be  similar  to 
the  TEG  pyrolysis  that  is  described  by  a  two  step  decomposition  process  using  two  generalized 
reaction  parameters  §f5  and  §f6  with  a  periodically  supplied  molar  concentration,  nTMj,  of  TMI. 
The  formation  of  GaP  and  InP  and  its  incorporation  in  the  underlying  film  is  summarized  in  the 
reactions  Eqs.  (9)  and  (10).  The  composition,  x,  for  the  compound  semiconductor  Gaj.xInxP  is 
expressed  as  the  averaged  ratio  of  molar  concentration  over  a  cycle  sequence 


and  the  instant  film  growth  rate  gp  is  given  by 


Sfl  - 


1 

A 


+ 


(12) 


At  this  point,  the  SRL  is  treated  as  a  homogeneous  ideal  solution  and  the  surface  area,  A,  is 
assumed  to  be  constant  for  simplicity.  Also  note  that  the  surface  structure,  number  of  reaction 
sides,  and  inhomogeneous  reactions  are  not  explicitly  addressed  at  this  point  and  are  integrated  into 
the  reaction  parameters  §f4  and  §f7 . 

The  temporal  thickness  evolution  of  the  SRL  is  given  by 


dl(t) 


+  n2  V2  +  n3  V3  +  n5  V5  + 


(13) 


where  Vj  are  the  molar  volumes  of  the  constituents  in  the  SRL. 

Based  on  the  above  reaction  chemistry,  we  model  the  linkage  between  the  measured  PR  signals 
and  the  surface  kinetics  on  the  basis  of  a  reduced  order  surface  kinetics  model  and  a  four-media 
stack:  ambient  /  SRL  /  epilayer  /  substrate,  which  represents  the  simplest  possible  description  of  the 
optical  response  under  the  conditions  of  PCBE  processes.  For  the  interpretation  of  the  time- 
dependence  of  the  four  media  stack  reflectance,  R4  (t),  in  terms  of  the  chemical  kinetics  in  the 
SRL  that  drives  epitaxial  growth,  the  dielectric  function  of  the  SRL,  8j,  must  be  linked  to  its 
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composition.  Such  a  linkage  can  be  established  by  approximating  the  dielectric  function  of  the  SRL 
through  an  effective  dielectric  function  ei,  parameterized  and  expressed  as  the  sum  over  all  molar 
fractions  Xj  contributing  to  the  SRL: 

8i(C0)  =  8oo+  Xxi(t)Fi(«)  andxj(t)=  .  (14) 

i#4,7  2j  knk(t) 

The  dielectric  function  of  the  SRL  is  obtained  by  summing  over  the  contributions  of  all  its 
constituents,  identified  by  the  label  i.  In  equation  (14),  F;(co)  denotes  optical  response  factors 

associated  with  transitions  that  characterize  specific  molecular  fragments,  and  co  denotes  the 

frequency  at  which  £1  is  evaluated. 

This  ROSK  model  provides  a  description  how  to  relate  changes  in  composition  and  thickness  of 
the  SRL  to  an  effective  dielectric  function  8i(co,t)  and  dj(t),  respectively.  It  allows  also  to  obtain 
the  instantaneous  composition,  x,  and  growth  rate  gfi(t)  of  the  GaxIii]_xP  film.  The  ROSK  data 
are  incorporated  in  Fresnel's  equation  that  determines  the  reflectance  amplitude,  rr,  of  the  p- 
polarized  light  as  follows.  Consider  the  four  layer  media  composed  of  ambient  /  SRL  /  film  / 
substrate.  We  model  the  reflection/refraction  of  the  surface  reaction  layer  by  an  effective  medium 
with  the  homogeneous  dielectric  function  8 j  (t)  and  the  thickness  dj(t).  Let  us  denote  the  four 
media  by  the  indices  n=0,l,2,3  labeled  from  the  ambient  to  the  substrate.  The  reflection  coefficient 
rn_i  n  from  the  (n-l)-st  media  to  n-th  media  is  given  by 


£n  V 

8n-t  -  £q  sin"(p  -  8n_!  -y 

.  ? 

£n  -  £q  sin  cp 

£n  A/ 

8 n-i  -  e0  sin2(p  +  8n_!  A 

/  8n  -  £0  sin2cp 

(15) 


where  8n  is  the  complex  dielectric  function  of  the  n-th  media.  The  factor  On  for  the  n-th  media  is 
given  by 

_  2  7t  dn  ,  re 

On  =  V  £n  -  e0  Sin  (p  , 

A.  (16) 

where  dn  is  the  thickness  of  the  n-th  media.  For  the  formulation  of  the  control  problem  we  analyze 
the  multi-layer  film  stack  of  GaxInj_xP  with  different  composition  x  using  the  virtual  interface  (V-I) 

method  described  by  Aspnes27-29  Consider  the  above  formulated  four  layer  media  composed  of 
ambient  /  SRL  /  film  /  substrate  and  replace  it  through  a  four  media  stack  that  is  built  up  by 

•  the  ambient  (0), 

•  the  surface  reaction  layer,  SRL,  represented  by  (1) 
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•  a  the  near  surface  layer  (u)  describing  the  growing  film  with  composition  x,  and 

•  a  virtual  substrate  (v),  represented  by  its  complex  virtual  reflectance  coefficients  rv. 
With  this,  the  total  reflectance  amplitude  is  given  as 


rr4 


r01  + 


1  +  r0j  f  e 


i/p-2i®i 


with  l / : 


-2  i 

rlu  +  rv  e 

1  ,  -2  i  O,, 

1  +  rlu  rv  e 


(17) 


Here,  iq | ,  r]u  and  are  functions  in  £j(t)  and  dj(t).  The  virtual  reflectance  coefficient  rv  can  be 
determined  by  knowing  the  thicknesses  and  dielectric  constants  for  each  of  the  underlying  layers. 
However,  even  under  exact  knowledge  of  each  of  these  parameters,  the  recursive  formula  that 
determines  rv  can  be  numerically  unstable  and  non-robust. 


r - \ 


Virtual  Reactor 


v _ / 

Fig.  5:  Control  of  heteroepitaxial  GalxInxP  growth:  The  Compensator  Design  consists  of  three 
elements:(l)  ROSKM  described  by  f;  (2):  Filter  gains  Gj(t)  based  on  Nonlinear-filtering 
techniques  and  (3)  Feedback  law  K  based  on  Dynamical  programming 


At  this  point,  we  can  establish  a  mathematical  control  loop,  as  schematically  outlined  in  Fig.  5, 
for  the  control  of  the  deposition  process  with  a  control  of  composition  x  in  the  GaxIiij.xP  surface 

layer,  utilizing  the  real-time  measured  PR  signals  at  the  two  different  angles  of  incidence  cpj,  i  =1,2 
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hj(  £j(t),  d,(t)  ,gr(t),  x(t)  )  +  noise 


Ri(t)  = 


(18) 


where  the  output  function  hj  are  determined  by  Eqs.  (15) -(18).  The  nonlinear  filtering  algorithms 

applied  for  real-time  estimates30  yield  the  composition  x,  the  growth  rate  gr(t)  and  the  dielectric 
function  of  the  surface  layer.  . 


Now  we  can  formulate  the  problem  of  controlling  the  growth  rate  and  composition  of  GaxIri|_xP 
as 


(19) 


subject  to  — n(t)  =  f(n(t))  +  B  s(t)  +  noise, 
dt 

where  gt[  is  the  desired  growth  rate  per  cycle  and  xt|  the  desired  composition. 


The  application  of  the  nonlinear  filtering  algorithm  is  demonstrated  in  Fig.  6  during  growth  of  a 

o 

parabolic  graded  500A  thick  GalnP  structure.  The  PR  responses  PR75  and  PR70  are  the  feedback 
signals  to  control  TEG  and  TMI  flows  based  on  Nonlinear-filtering  techniques  and  Dynamical 
programming  shown  in  Fig.  5.  In  this  experiment,  the  filtering  output  was  integrated  over  10 

o 

cycle  sequences  (60  sec,  20A  resolution)  to  adjust  the  control  signals  for  the  TEG  and  TMI  flows. 
The  high  integration  time  caused  some  controller  instabilities,  observed  in  the  setting  of  the 
TMLTEG  flow  ratio  (see  fluctuations  in  flow  ratio  control  signal  in  Fig.  6),  which  shows  a 
stepwise  adjustment  /  fluctuation  in  the  TEG  and  TMI  flows.  The  two  insets  in  Fig.  6  show 
enlarged  the  evolutions  of  the  PR  signals  (30  sec),  taken  for  the  same  time.  The  different 
amplitudes  are  related  to  the  distinct  composition  and  thickness  of  the  SRL  as  described  by  the 
ROSK  model. 


Figure  7  shows  a  offline  analysis  and  parameter  estimation  of  the  PR  signals  during  growth  of 
GaP  on  Si(001)  at  350  °C.  The  experimental  PR  signals  were  recorded  for  two  angles  of 
incidence,  tp=71°  andtp=75.2°,  at  A=650±5nm  and  at  k=632.8  nm,  respectively.  The  two  insets  in 

Fig.  7  show  enlarged  the  evolutions  of  the  PR  signals  at  different  times  during  growth  of  GaP  on 
Si(001)  and  are  compared  with  the  simulation.  The  off-line  simulated  responses  using  the  ROSK 
model  and  nonlinear  filtering  are  shown  with  an  offset  in  order  to  compare  it  with  the  experiment. 
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time  (s) 

o 

Fig.  6:  Closed-loop  feedback  controlled  growth  of  a  500A  parabolic  Gal  xInxP  heterostructure. 


The  steps  of  the  generation  of  a  set  of  simulated  data  are  shown  in  detail  in  Fig.  8,  for  the  TEG 
pulse  1.3  -  1.6  s  (and  TBP  pulse  0.0  -  0.8  s)  data.  The  three  SRL  components  are  the  result  of  the 
source  pulses  and  the  ROSK  reaction  approximations.  From  the  SRL  components,  the  SRL 
thickness  and  dielectric  function  are  found.  These  values  then  contribute  to  the  calculated 
reflectance.  The  fit  of  the  simulated  fine  structure  to  the  experimental  data  as  shown  in  Fig.  8(e)31. 
The  parameters  found  by  the  minimization  resulted  in  efilm=10.6  -  0.06i;  £Sllbstrate  =  15.82  -  0.27i;  an 
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average  SRL  dielectric  function  of  16.82-4.47/  ;  and  an  average  film  growth  rate  of 
0.365 angstroms  I  s  . 


Figure  7:  Online  state  estimation  of  surface  kinetics  based  on  two  PRS-signals. 


A  change  in  cycle  sequence,  flow  ratio,  flux,  temperature  or  total  system  pressure,  leads  to 
different  results,  with  no  predictable  correlation  found  yet.  Presently,  most  of  the  relevant 
reactions  pathways  and  their  rate  constants  are  not  known,  with  almost  no  information  on  reactions 
kinetic  data  at  growth  temperature.  Accessing  these  data  will  require  real-time  infra-red  PRS  to 
analyze  molecular  specific  absorption  signatures  and  their  correlation  with  real-time  mass 
spectroscopic  signals.  Such  an  analysis  would  allow  a  gain  a  better  understanding  of  the  surface 
reactions  kinetics,  the  validation  of  the  ROSK  model,  and  would  provide  quantitative  input 
parameter  for  more  detailed  modeling  of  the  growth  process. 
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Figure  8: 


TBP 


Simulation  of  surface  reaction 
kinetics  and  PR  response 
during  heteroepitaxial  GaP 
growth  on  Si  under  pulsed 
organometallic  precursors 
exposure  of  TBP  and  TEG. 
The  molar  concentrations  of 
the  surface  constituents,  their 
effective  surface  layer 
thickness  and  the  dielectric 
optical  properties  are 
calculated  using  equations  (5) 
through  (15). 

The  simulated  PR  response  is 
compared  with  experimental 

results  obtained  at  cp=75.1°, 
X=632.8  nm. 


Summary 

We  reported  the  compositionally  and  thickness  controlled  growth  of  GalnP  heterostructures, 
using  PRS  as  feedback  sensor.  An  on-line  parameter  estimate  of  the  state  functions  iV(t)  using  the 
PRS  signals  is  used  to  analyze  the  PRS  fine  structure  and  to  provide  the  control  signals  s(t)  for 
TEG  and  TMI  flows.  We  introduced  a  reduced  order  kinetics  (ROSK)  model  using  generalized 
reaction  rate  parameters  to  describe  the  thickness  and  composition  evolution  of  the  SRL.  The  molar 
concentrations  nj  in  the  SRL  are  directly  linked  to  instantaneous  growth  rate  and  composition  of  the 
underlying  film. 
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